---------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  /Users/wesleyblundell/Dropbox/Wes_Supplemental_NEI/JAERE_Replication/Replication_package/log/Table_A5_Results.log
  log type:  text
 opened on:  15 May 2024, 11:40:42

. eststo clear

. use "$rootpath/processed_data/Estimation_Data.dta", clear

. 
. reghdfe  log_damage55 treatXpost if year>2009, absorb(facility_id n5#year) vce(cluster facility_id)
(dropped 2885 singleton observations)
(MWFE estimator converged in 250 iterations)

HDFE Linear regression                            Number of obs   =    813,477
Absorbing 2 HDFE groups                           F(   1, 108733) =      34.02
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9468
                                                  Adj R-squared   =     0.9380
                                                  Within R-sq.    =     0.0002
Number of clusters (facility_id) =    108,734     Root MSE        =     0.3593

                      (Std. Err. adjusted for 108,734 clusters in facility_id)
------------------------------------------------------------------------------
             |               Robust
log_damage55 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  treatXpost |   -.057528   .0098628    -5.83   0.000     -.076859   -.0381971
       _cons |    .829215   .0016336   507.59   0.000     .8260132    .8324169
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
 facility_id |    108734      108734           0    *|
     n5#year |      7170           1        7169     |
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

. estimates store outlier_damage

. 
. //Emissions results
. reghdfe  log_emissions4 treatXpost if year>2009, absorb(facility_id n5#year) vce(cluster facility_id)
(dropped 2885 singleton observations)
(MWFE estimator converged in 254 iterations)

HDFE Linear regression                            Number of obs   =    813,477
Absorbing 2 HDFE groups                           F(   1, 108733) =      38.46
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9477
                                                  Adj R-squared   =     0.9390
                                                  Within R-sq.    =     0.0002
Number of clusters (facility_id) =    108,734     Root MSE        =     0.3074

                      (Std. Err. adjusted for 108,734 clusters in facility_id)
------------------------------------------------------------------------------
             |               Robust
log_emissi~4 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  treatXpost |  -.0547269   .0088243    -6.20   0.000    -.0720225   -.0374314
       _cons |   .6927449   .0014616   473.96   0.000     .6898801    .6956096
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
 facility_id |    108734      108734           0    *|
     n5#year |      7170           1        7169     |
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

. estimates store outlier_emission

. 
. drop if key_estimation_sample==0
(273,947 observations deleted)

. sort facility_id year

. egen num_obs=nvals(year), by(facility_id)

. 
. //Balanced sample damages. 
. reghdfe  log_damage55 treatXpost if num_obs==10, absorb(facility_id n5#year) vce(cluster facility_id)
(dropped 550 singleton observations)
(MWFE estimator converged in 183 iterations)

HDFE Linear regression                            Number of obs   =    336,280
Absorbing 2 HDFE groups                           F(   1,  33668) =      24.28
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9408
                                                  Adj R-squared   =     0.9330
                                                  Within R-sq.    =     0.0004
Number of clusters (facility_id) =     33,669     Root MSE        =     0.4230

                       (Std. Err. adjusted for 33,669 clusters in facility_id)
------------------------------------------------------------------------------
             |               Robust
log_damage55 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  treatXpost |  -.1122007   .0227684    -4.93   0.000    -.1568275   -.0675739
       _cons |   1.319793   .0040705   324.24   0.000     1.311814    1.327771
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
 facility_id |     33669       33669           0    *|
     n5#year |      5655           1        5654     |
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

. estimates store balanced_damage

. 
. ///This is balanced sample emissions
> reghdfe  log_emissions4 treatXpost if num_obs==10, absorb(facility_id n5#year) vce(cluster facility_id)
(dropped 550 singleton observations)
(MWFE estimator converged in 183 iterations)

HDFE Linear regression                            Number of obs   =    336,280
Absorbing 2 HDFE groups                           F(   1,  33668) =      32.84
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9440
                                                  Adj R-squared   =     0.9365
                                                  Within R-sq.    =     0.0006
Number of clusters (facility_id) =     33,669     Root MSE        =     0.3545

                       (Std. Err. adjusted for 33,669 clusters in facility_id)
------------------------------------------------------------------------------
             |               Robust
log_emissi~4 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  treatXpost |  -.1142533   .0199381    -5.73   0.000    -.1533327   -.0751739
       _cons |   1.072635   .0035645   300.92   0.000     1.065649    1.079622
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
 facility_id |     33669       33669           0    *|
     n5#year |      5655           1        5654     |
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

. estimates store balanced_emission

. 
. //This is extreme damages outliers dropped, 100 million
. reghdfe  log_damage55 treatXpost if key_estimation_sample==1&tot_damage<1e8, absorb(facility_id n5#year) vce(cluster facility_id)
(dropped 19 singleton observations)
(MWFE estimator converged in 240 iterations)

HDFE Linear regression                            Number of obs   =    697,655
Absorbing 2 HDFE groups                           F(   1,  81363) =      49.74
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9392
                                                  Adj R-squared   =     0.9304
                                                  Within R-sq.    =     0.0003
Number of clusters (facility_id) =     81,364     Root MSE        =     0.3582

                       (Std. Err. adjusted for 81,364 clusters in facility_id)
------------------------------------------------------------------------------
             |               Robust
log_damage55 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  treatXpost |  -.0695985   .0098682    -7.05   0.000    -.0889401   -.0502569
       _cons |   .8551122   .0016912   505.64   0.000     .8517975    .8584268
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
 facility_id |     81364       81364           0    *|
     n5#year |      6831           1        6830     |
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

. estimates store no_outlier_damage

. 
. //100 million, emissions
. reghdfe  log_emissions4 treatXpost if key_estimation_sample==1&tot_damage<1e8, absorb(facility_id n5#year) vce(cluster facility_id)
(dropped 19 singleton observations)
(MWFE estimator converged in 242 iterations)

HDFE Linear regression                            Number of obs   =    697,655
Absorbing 2 HDFE groups                           F(   1,  81363) =      59.63
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9407
                                                  Adj R-squared   =     0.9321
                                                  Within R-sq.    =     0.0004
Number of clusters (facility_id) =     81,364     Root MSE        =     0.3048

                       (Std. Err. adjusted for 81,364 clusters in facility_id)
------------------------------------------------------------------------------
             |               Robust
log_emissi~4 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  treatXpost |  -.0675649   .0087496    -7.72   0.000    -.0847142   -.0504157
       _cons |   .7048204   .0014995   470.05   0.000     .7018815    .7077594
------------------------------------------------------------------------------

Absorbed degrees of freedom:
-----------------------------------------------------+
 Absorbed FE | Categories  - Redundant  = Num. Coefs |
-------------+---------------------------------------|
 facility_id |     81364       81364           0    *|
     n5#year |      6831           1        6830     |
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

. estimates store no_outlier_emission

. 
. //Main results six digit naics. 
. reghdfe  log_damage55 treatXpost if key_estimation_sample==1, absorb(facility_id naics#year) vce(cluster facility_id)
(dropped 976 singleton observations)
(MWFE estimator converged in 230 iterations)

HDFE Linear regression                            Number of obs   =    698,037
Absorbing 2 HDFE groups                           F(   1,  81366) =      39.16
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9418
                                                  Adj R-squared   =     0.9330
                                                  Within R-sq.    =     0.0002
Number of clusters (facility_id) =     81,367     Root MSE        =     0.3600

                       (Std. Err. adjusted for 81,367 clusters in facility_id)
------------------------------------------------------------------------------
             |               Robust
log_damage55 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  treatXpost |  -.0625375   .0099941    -6.26   0.000     -.082126    -.042949
       _cons |   .8673023   .0017128   506.36   0.000     .8639452    .8706594
------------------------------------------------------------------------------

Absorbed degrees of freedom:
------------------------------------------------------+
  Absorbed FE | Categories  - Redundant  = Num. Coefs |
--------------+---------------------------------------|
  facility_id |     81367       81367           0    *|
   naics#year |      9907           1        9906     |
------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

. estimates store naics6_damage

. 
. //Main results six digit naics. 
. reghdfe  log_emissions4 treatXpost if key_estimation_sample==1, absorb(facility_id naics#year) vce(cluster facility_id)
(dropped 976 singleton observations)
(MWFE estimator converged in 232 iterations)

HDFE Linear regression                            Number of obs   =    698,037
Absorbing 2 HDFE groups                           F(   1,  81366) =      50.57
Statistics robust to heteroskedasticity           Prob > F        =     0.0000
                                                  R-squared       =     0.9429
                                                  Adj R-squared   =     0.9343
                                                  Within R-sq.    =     0.0003
Number of clusters (facility_id) =     81,367     Root MSE        =     0.3051

                       (Std. Err. adjusted for 81,367 clusters in facility_id)
------------------------------------------------------------------------------
             |               Robust
log_emissi~4 |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
  treatXpost |  -.0627157   .0088193    -7.11   0.000    -.0800016   -.0454299
       _cons |   .7138576   .0015115   472.29   0.000     .7108951      .71682
------------------------------------------------------------------------------

Absorbed degrees of freedom:
------------------------------------------------------+
  Absorbed FE | Categories  - Redundant  = Num. Coefs |
--------------+---------------------------------------|
  facility_id |     81367       81367           0    *|
   naics#year |      9907           1        9906     |
------------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation

. estimates store naics6_emission

. 
. 
. esttab naics6_damage outlier_damage balanced_damage no_outlier_damage naics6_emission outlier_emission balanced_emission no_outlier_emission using "$ro
> otpath/table_graph/TBL_A5.tex",  star( * 0.1 ** 0.05 *** 0.01) cells(b(star fmt(3)) se(par fmt(3))) keep(treatXpost _cons) scalars(N F) ar2 replace 
(note: file /Users/wesleyblundell/Dropbox/Wes_Supplemental_NEI/JAERE_Replication/Replication_package/table_graph/TBL_A5.tex not found)
(output written to /Users/wesleyblundell/Dropbox/Wes_Supplemental_NEI/JAERE_Replication/Replication_package/table_graph/TBL_A5.tex)

. 
. log close
      name:  <unnamed>
       log:  /Users/wesleyblundell/Dropbox/Wes_Supplemental_NEI/JAERE_Replication/Replication_package/log/Table_A5_Results.log
  log type:  text
 closed on:  15 May 2024, 11:44:14
---------------------------------------------------------------------------------------------------------------------------------------------------------
